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We report transport studies on the 5 nm thick Bi 2 Se 3 topological insulator films which are grown via molecular 
beam epitaxy technique. The angle-resolved photoemission spectroscopy data show that the Fermi level of the 
system lies in the bulk conduction band above the Dirac point, suggesting important contribution of bulk 
states to the transport results. In particular, the crossover from weak antilocalization to weak localization in the 
bulk states is observed in the parallel magnetic field measurements up to 50 Tesla. The measured 
magneto-resistance exhibits interesting anisotropy with respect to the orientation of parallel magnetic field B 7/ 
and the current I, signifying intrinsic spin-orbit coupling in the Bi 2 Se 3 films. Our work directly shows the 
crossover of quantum interference effect in the bulk states from weak antilocalization to weak localization. It 
presents an important step toward a better understanding of the existing three-dimensional topological 
insulators and the potential applications of nano-scale topological insulator devices. 

Topological insulators (TIs), a new type of quantum materials characterized by the existence of a gap in the 
energy spectrum and the symmetry-protected surface states within the bulk gap, have attracted significant 
interest in recent years 1 " 3 . As a result of the strong spin-orbit coupling (SOC), the spin-momentum locked 
surface states always show weak antilocalization (WAL) effect 4 " 6 due to the acquired n Berry phase when 
completing a time-reversed self-crossing path. Plentiful interesting behaviors of TIs, particularly in the topo- 
logically protected surface states, have been explored through electrical transport experiments 7 " 22 . Whereas, it is 
known that the bulk states of the existing three-dimensional (3D) TIs are not perfectly insulating even at low 
temperatures, and as a direct consequence, the bulk states can also influence the transport behaviors of TIs. In 
order to develop a comprehensive understanding on the exotic properties of TIs, it is important to systematically 
study the transport of TIs taking into account the effects from the bulk states. 

Contrasting to the symmetry protected surface states, the TI bulk states with strong SOC can lead to either 
WAL or weak localization (WL) depending on the experimental conditions 5,6 ' 23 . Recently, the WL effect of bulk 
states has been studied in the ultrathin Bi 2 Se 3 films under perpendicular magnetic field 24 . However, it is difficult to 
distinguish the bulk state features from the surface state behaviors by the perpendicular field transport measure- 
ments. In comparison, the parallel field magneto -resistance (MR) properties of TI films reveal more bulk-state 
information 25 ' 26 , and therefore, the parallel magnetic field transport measurements provide a powerful tool in 
investigating the quantum interference effects in the TI bulk states. 

Here we comprehensively study the transport properties of 5 quintuple layers (QLs) Bi 2 Se 3 films grown on epi- 
ready a-Al 2 0 3 (0001) substrates by molecular beam epitaxy (MBE). The in situ angle-resolved photoemission 
spectroscopy (ARPES) result shows that two quantum well states of the conduction band, as well as the surface 
states, contribute to the density of states near the Fermi level. The MR is measured in pulsed high magnetic field 
(PHMF) up to 50 Tesla (T). When the external field is perpendicular to the film plane, linear MR is observed 
above 16 T. In parallel magnetic field, MR behaviors of the TI films exhibit interesting anisotropy when the field 
orientation is tuned between B /7 J_ I and B /7 // 1 configurations. Specifically, while the field is increased up to 50 T, 
the MR always exhibits a distinctive switch from positive to negative when B// _L I; whereas, when B// // 1, such 
switch either occurs at a much larger field, or does not occur at all within 50 T. The theoretical analysis and the 
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quantitative fitting for the experimental results indicate that such 
parallel field MR behaviors of the ultrathin TI films can be well 
explained using the WAL-WL crossover mechanism in the TI bulk 
states. To be specific, as a consequence of the strong SOC effect in the 
TI films, the phase coherence time generally presents a larger time 
scale than the intrinsic spin-flip time x S o- As a result, the WAL effect 
always occurs when the magnetic field is relatively weak, whereas, the 
WL effect arises when the field becomes sufficiently strong. From this 
perspective, the MR switch from positive to negative can be under- 
stood as a manifestation of the WAL-WL crossover in the TI bulk 
states. In addition, we observe novel MR anisotropy with respect to 
the relative orientation of the parallel magnetic field and the current, 
which can be qualitatively explained using the SOC mechanism. 
Combining the ARPES data, the transport results and the theoretical 
analysis, we reveal the crossover of quantum interference effect, from 
WAL to WL, in the bulk states of the ultrathin Bi 2 Se 3 films. 

Results 

The 5 QLs Bi 2 Se 3 films are epitaxially grown on sapphire (0001) 
substrates in an ultra-high vacuum MBE-ARPES-STM (scanning 
tunneling microscope) combined system 27 . Figure 1(a) shows the 
reflection high- energy electron diffraction (RHEED) pattern of a 
typical MBE-grown 5 QLs Bi 2 Se 3 film, where the sharp 1X1 streaks 
demonstrate that the film has good single crystal quality. The band 
map of the film taken in T-K direction is shown in Figure 1(b). In the 
ARPES measurement, photoelectrons are excited by an unpolarized 
He-I light (21.21 eV), and collected by a Scienta SES-2002 analyzer. 
The topological surface states are observed within the bulk gap. A 
quite small gap is opened in the Dirac cone just below the Dirac point 
due to the coupling between the top and bottom surface states. 
Moreover, as shown in Figures 1(b) and (c), the quantum well states 
of the conduction band have a significant contribution to the Fermi 
surface. These observations indicate that our films are TIs with the 
Dirac surface states and the bulk states play a great role in the charge 
transport. 

In order to explore the transport properties of the films by ex situ 
measurements, 20 -nm- thick insulating amorphous Se is deposited 
on the TI films as a protection capping layer (Figure 1(d)). The Hall 
bars with channel dimensions of 600 X 400 um 2 are fabricated using 
standard photo -lithography method. An optical image of the typical 
Hall bar structure is shown in the inset of Figure 2(a). The transport 



properties of the films were measured mainly by using the standard 
four-electrode method in a 16 T Quantum Design PPMS (physical 
property measurement system) and PHMF at the Wuhan National 
High Magnetic Field Center. The non- destructive pulsed magnet has 
a rise time of 15 ms and fall time of 135 ms. The field and current 
directions are reversed in the measurements to improve the precision 
of the results by eliminating the parasitic voltages. Typical results 
from the sample 1 (a 5 QLs Bi 2 Se 3 film) are shown in the main text. 
(The results are repeatable in a collection of samples.) 

The dependence of the sheet longitudinal resistance R 0 on the 
temperature was measured firstly (for the sample 1), which is shown 
in Figure 2(a). The resistance decreases almost linearly when the 
temperature decreases from 300 K to about 50 K. Whereas, when 
the temperature drops approximately below 12 K, the resistance 
develops an upturn. The upturn resistance, as shown in 
Figure 2(b), exhibits a logarithmic dependence on the temperature, 
which can be understood in terms of the quantum correction 
induced by the electron- electron interaction 13 . The resistance upturn 
becomes larger at 8 T than that at 0 T, which is due to the suppres- 
sion of the WAL effect in the presence of the perpendicular magnetic 
field resulting in enhanced resistance. The Hall trace of the sample 1 
at 4.2 K is shown in Figure 2(c). Utilizing the slope at low fields, we 
calculate the sheet carrier density n s ~ 1.89 X 10 13 cm" 2 and mobil- 
ity \i ~ 316.6 cm 2 /(V*s) at 4.2 K. 

Next, the MR behavior of the sample 1 at relatively low perpen- 
dicular fields (<16 T) was measured in PPMS and PHMF, respect- 
ively. The results in both measurements are consistent, which are 
shown in Figure 2(d). Specifically, these MR curves show typical dips 
around 0 T at low temperatures as a result of the WAL effect in 
TIs 7,12 " 15 , which can be clearly seen in the inset of Figure 2(d). In 
addition, a linear MR appears above 16 T at both 4.2 K and 77 K. 
This is consistent with the previous observation in a 5 QLs Bi 2 Se 3 film 
grown on the same substrate by MBE 15 . Note that similar linear 
behavior of MR has been observed in many TI systems 28 " 35 , yet, the 
underlying origin remains to be clarified. Some attribute such linear 
MR to a quantum origin, including the Abrikosov's quantum linear 
MR model which relates the linear MR to the linear energy dispersion 
of the gapless topological surface states 28 " 30 , and the Hikami-Larkin- 
Nagaoka theory which interprets the linear-like MR at high fields in 
terms of the WAL mechanism 31,32 . Other theories associating the 
linear MR with a classical origin have also been proposed, such as 




Figure 1 | In situ characterizations of the 5QLs Bi 2 Se 3 thin film grown by MBE and schematic structure for ex situ transport measurement, (a) RHEED 
pattern with sharp 1X1 streaks, (b) ARPES band spectra along the T-K direction. The white dashed line indicates the Fermi level. The arrows label the 
quantum well states, the Dirac point and the small gap, respectively, (c) The corresponding momentum distribution curves of (b). The blue line 
indicates the Fermi level, (d) The schematic structure for transport measurements. The thickness is not to scale. 
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Figure 2 | Transport properties of sample 1. (a) The sheet resistance R 0 versus temperature with an upturn in low temperature regime. The inset is 
an optical image of the Hall bar structure, (b) Normalized upturn resistances at fixed perpendicular fields, (c) Hall trace at 4.2 K in a magnetic field up to 
50 T. (d) Normalized perpendicular field MR. The inset shows the MR dips around 0 T which are attributed to the WAL effect at low temperatures. 
The broken lines are guides to the eyes. 



the inhomogeneity model by Parish and Littlewood which stresses 
the influence from the sample inhomogeneity 34 ' 35 . 

Interestingly, a different set of rich MR behavior is observed when 
a magnetic field parallel to the film plane is applied. In particular, as 
clearly shown in Figure 3, the normalized MR of sample 1 along the 
[110] direction in parallel field, exhibits a significant dependence on 
the relative orientation of B// and I. The behavior of MR when the 
parallel field is oriented along B// J_ I (the field increases up to 50 T) 
is shown in Figure 3(a). At 4.2 K, the MR first increases to a max- 
imum at about 25 T, and then decreases as a downward parabolic- 
like curve. When the temperature is increased to 77 K, the MR 
slightly increases within 1 1 T and then decreases monotonically up 
to 50 T. In comparison, the MR behavior for the field orientation B// 
// 1 (Figure 3(b)) is different, only showing positive MR within 50 T 
at both 4.2 K and 77 K. 

To further understand the parallel field MR behavior, we also 
measured the parallel field MR along the [ — 110] direction of the 
sample 1. The results are shown in Figures 4(a) and (b). The quasi- 



four-electrode method (insets of Figures 4(a) and (b)) is used only 
here due to the restriction of sample size. It is worth mentioning that 
such measurement method can introduce contact resistance that can 
lead to a weaker MR response without changing the curve shape 
(Supplementary Information). In this case, the MR switch from pos- 
itive to negative is also clearly observed at 4.2 K and 77 KwhenB// J_ 
I, but much weaker when B// // 1. For a fixed temperature, the switch- 
ing field for B// // 1 is larger than that for B f/ J_ I. In order to explore 
the general dependence of the MR switch on the B /r I orientation, the 
parallel field MR of sample 1 with the current along the [110] dir- 
ection under different position 0 was measured (Figure 4(c)). The 
sample position 0 refers to the angle between B// and the normal 
direction of I along the film plane (the inset of Figure 4(c)). Here 0 = 
0° means B /7 J_ I and 0 = 90° means B /7 // 1. When the magnetic field 
increases up to 16 T, the MR at 100 K increases when B /7 // I, but 
decreases when B /7 J_ I. When 0 = 45°, the MR seems to show a 
competition between positive and negative MR. Normalized MR at 
100 K as a function of 0 when B// = 15 T is shown in Figure 4(d). 




B,(T) B,(T) 

Figure 3 | Normalized MR of sample 1 in parallel field along [110] direction, (a) B// _L I and (b) B// // 1. The triangles are experimental data and 
the black dashed lines are fitting curves by Equation ( 1 ) . The MR behaviors can be well explained by the WAL- WL crossover mechanism in TI bulk states. 
The insets show the corresponding measurement configurations. [110] and [ — 110] are two different crystal directions along the film plane. 
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Discussion 

A positive-to-negative MR switch of the ultrathin Bi 2 Se 3 films is 
observed when the parallel magnetic field increases. The switch 
occurs in a broad range of temperature and exhibits anisotropy when 
the relative orientation of B /f and I is varied. The parallel field MR 
data exhibit tiny variation in the presence of a large magnetic field, 
indicating a quantum correction regime. Taking into account the fact 
that the Fermi level locates at the bulk conduction band, and that the 
film shows quite large carrier density and lower mobility compared 
to those reported for the surface states of TI 10 , we conclude that the 
major contribution to the charge transport in our films comes from 
the bulk states instead of the surface states. We emphasize that, while 
similar negative MR behavior observed in Sn-doping Bi 2 Te 3 films 26 
under parallel field was reported, a clear interpretation of the MR 
data and the comparison of different behaviors in Sn-doping samples 
are still missing. In contrast, for our present experimental results, the 
parallel field MR switch phenomenon of Bi 2 Se 3 films can be well 
explained in terms of the quantum interference effects in the TI bulk 
states. In particular, as we shall describe in detail below, the positive 
and negative MR can be attributed to the presence of WAL and WL, 
respectively. 

For a TI system with strong SOC, the correction to the MR beha- 
vior due to the effect of quantum interference 6 is characterized by two 
time scales: the dephasing time and the spin-flip time x so . In the 
regime > x so , the spin-orbit scattering is significant, which gives 
rise to frequent spin flips. As a result, the quantum interference 
between the time reversed trajectories is destructive, leading to the 
occurrence of WAL and a positive MR. In the regime where x^ is 
comparable to x S o or x 9 < x S o> the spin-orbit scattering effect is weak 
and the scattering process barely affects the spin orientation. In such 
case, negative MR occurs with increasing magnetic field as a WL 
feature. In the intermediate regime where x^ > x so > MR first 
increases to a peak value and then decreases, signifying a crossover 
from WAL to WL. Overall, the quantum interference behavior of a 
system is crucially dependent on the ratio between x 9 and x so . 



In our thin films, the mean free path / can be estimated according 
to / = kkpfi/e, where k F can be obtained through the sheet density n s . 
For the sample 1 described here, / is approximately 20 nm, much 
larger than the thickness d = 5 nm. The parallel field transport 
thereby falls into the Dugaev-Khmelnitskii (DK) regime where the 
resistance change is given by 36,37 



— = -N 
R 1 2n 2 h 



3 / L 2 



In 1 + 



(1) 



Here AR = R(B) — R(B = 0). R can be regarded as the resistance in 
the absence of external magnetic field due to the quite small vari- 
ation, and R D is the sheet resistance. N A is the number of channels 
contributing to the charge transport. In addition, three length scales 
are involved: L so is the spin-flip length (the distance travelled by an 
electron before its spin direction is changed by the scattering, which 
depends on the strength of SOC in the bulk states. Please see 
Supplementary Information.), L^ is the phase coherence length, 

and L|| =Al} B ^/TJd d with L B being the magnetic length 



L B = \JhJeB. By varying only two fitting parameters L 9 and L so , 
we present a quantitative fitting of the experimental data in 
Figure 3. In the case of B// _L I, we have L^ ~ 140 nm, L so ~ 
30 nm at 4.2 K and L^ ~ 40 nm, L so ~ 20 nm at 77 K. And for 
B// // 1, L 9 ~ 140 nm, L so — 18 nm at 4.2 K. Based on these fitting 
results, L 9 ~ 40 nm and L so —14 nm at 77 K when B// // 1 can be 
obtained by calculation (see Supplementary Information for details). 
Substituting these calculated values into Equation (1), we quantita- 
tively reproduce the MR result at 77 K in B// // I configuration as 
shown in Figure 3(b). The tiny mismatches of the fitting MR curves 
may arise from the neglect of the electron- electron interaction 13 , 
which can actually play a role in the MR behavior. 

It is noted that N f = 2 in the fitting is identical to the number of 
bulk bands at the Fermi surface shown in Figure 2(b). While we 
cannot fully determine whether the two channels are completely 
from the two subbands of the bulk state or from the hybridized 
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Figure 4 | Control experiments with sample 1. Normalized MR along [ — 110] direction when (a) B /7 _L I and (b) B// // 1. Normalized MR along [110] 
direction at 100 K when (c) 0 = 0°, 45°, 90° and (d) 0 continuously varies from 0° to 90°. 0 = 0° means B„ _L I and 0 = 90° means B„ // 1. The insets show 
the corresponding measurement configurations. 
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channels of bulk state and surface state, it is evident that the bulk state 
plays dominant role in the emergence of crossover behavior reported 
here. All the values of parameters shown above are fitting results and 
reasonable. The phase-coherence length L 9 decreases with temper- 
ature and L 9 is much larger than the mean free path as reported 
earlier 15 . Besides, the frequent scattering due to the edges in the 
ultrathin films makes L so comparable to the mean free path. Based 
on the fitting results of experimental data, L so is affected by the 
relative orientation of the parallel magnetic field and the current. 

Through the parallel field MR curves of the ultrathin Bi 2 Se 3 films, a 
clear image of the quantum interference effects with the WAL-WL 
crossover in the TI bulk states is depicted. Three cases are presented 
for the data in Figure 3 and 4. (I) For the B /f J_ I configuration, the 
system usually falls into the intermediate regime (x^ > x so ) at rela- 
tively low temperatures. As a result, WAL effect characteristically 
arises when the field is weak, whereas, when the magnetic field 
becomes sufficiently large, a crossover into WL occurs. The data of 
Figure 3(a) and 4(a) belong to this case and the crossover is observed 
within 50 T. (II) For the B /f II I configuration and relatively low 
temperatures, compared to (I), a weaker effect of the WAL-WL 
crossover is observed at a larger switching field due to a larger ratio 
tqAso as shown in Figure 4(b). In Figure 3(b), only the WAL effect is 
observed at 4.2 K and 77 K. This may indicate that the switching 
fields from WAL to WL in both conditions are larger than 50 T. (Ill) 
For either of the two geometries, the MR may show only WL effect at 
high temperatures when decreases and becomes comparable to 
x S c» even x^ < x S o The 0 = 0° data at 100 K in Figure 4(c) fall into 
this situation. 

The observed anisotropy feature of MR can be related to many 
factors, including the band structure asymmetry, crystalline aniso- 
tropy, and the relative orientation of the parallel magnetic field and 
the current. However, generally speaking, both the band asymmetry 
and the crystal anisotropy could result in a MR behavior that varies 
along different crystal directions. Since the anisotropy of MR along 
the [ — 110] direction remains consistent with that of the [110] dir- 
ection, we consider the relative directions of the parallel magnetic 
field and the current as the leading factor that contributes to the novel 
MR anisotropy reported here. The underlying physics can then be 
qualitatively understood in terms of the SOC mechanism. Consider a 
moving electron, which experiences an effective magnetic field B e 
that is perpendicular to the electron momentum. When B// ± I, the 
main transport path of the electron is perpendicular to the external 
field, in which the spin-flip process due to the SOC effect is strongly 
suppressed by the external magnetic field, leading to an effectively 
larger x S o- In distinct contrast, when B// // 1, the spin-flip process is 
hardly affected. This anisotropy in the effective SOC effect provides a 
mechanism that leads to the observed anisotropic MR here. Such 
explanation is also supported by Figures 4(c) and (d), which show 
an intimate relation between the parallel field MR and the sample 
position 0. Particularly, the negative MR behavior when B// _L I is 
consistent with a larger effective x so , while the positive MR when B /7 
// 1 agrees with a smaller x so . For arbitrary 0, the effective x so is given 
by the combination of both effects, which leads to the MR-0 char- 
acteristics shown in Figure 4(d). 

Interestingly, the anisotropy of the parallel field MR in thicker TI 
films (e.g. 45 QLs and 200 QLs Bi 2 Se 3 ) 38 seems to be opposite to the 
experimental results of the 5 QLs Bi 2 Se 3 films. Additionally, there is 
one order of magnitude difference of the MR amplitudes for the 
thicker films 38 between B//J_I and B////I, while it is not the case for 
the results presented here. In light of many different factors between 
the thick and ultrathin films, such as the thickness and the weight of 
bulk states at the Fermi surface, we are not able to achieve a unified 
understanding. Further theoretical investigations in this context 
would be valuable towards a quantitative comprehension of the 
anisotropic parallel field MR of TI films in various field- current 
orientations. 



In summary, the crossover from WAL to WL in TI bulk states is 
revealed clearly by the parallel field MR behaviors of the 5 QLs 
ultrathin Bi 2 Se 3 films. The WAL-WL crossover is largely affected 
by the relative orientation of the parallel magnetic field and the 
current, which can be understood in terms of an orientation-depend- 
ent effective spin-flip time. To be specific, when B// _L I, the MR 
within 50 T always exhibits the crossover at 4.2 K and 77 K, whereas, 
the crossover for B /7 // I occurs at a much higher magnetic field 
(possibly above 50 T), due to the anisotropic SOC effect under strong 
magnetic field. The novel transport properties of the ultrathin TI 
films demonstrated in this work are quite significant for a better 
understanding of the existing 3D TI materials and for the future 
potential magnetic applications of nano-scale TI devices. 

Methods 

Sample preparation. The high quality 5 QLs Bi 2 Se 3 films studied here were grown on 
oc-Al 2 0 3 (sapphire) (0001) substrates in an ultra-high vacuum MBE-ARPES-STM 
combined system with a base pressure better than 2 X 10" 10 mbar. Before sample 
growth, the sapphire substrates were first outgassed at 650°C for 90 min and then 
heated at 850° C for 30 min. High-purity Bi (99.9999%) and Se (99.999%) were 
evaporated from standard Knudsen cells. To reduce Se vacancies in Bi 2 Se 3 , the growth 
was kept in Se-rich condition with the substrate temperature at ~220°C. For the 
transport measurement, 20 nm Se was capped on the TI film as protection layer. 

Device fabrication and measurement. Hall bar structure with channel dimensions of 
600 X 400 urn 2 was fabricated using standard photo-lithography method. Positive 
photoresist S1813 was spun at 4000 rpm for 45 seconds on Bi 2 Se 3 film, followed by 
1 10°C baking for 60 seconds. With a mask of Hall bar pattern, the photoresist-coated 
sample was exposed to ultraviolet light (365 nm wavelength) with exposure power of 
8 mW/cm 2 for 7 seconds. The exposed part of photoresist was removed after 40 
seconds of developing (MFCD-26 developer). Then, bared area of Bi 2 Se 3 film with no 
photoresist on was wet-etched with 1 g of potassium dichromate in 10 ml of sulfuric 
acid and 20 ml of DI water. The expected etch rate for Bi 2 Se 3 was about 60 nm per 
minute, and DI water rinsing was needed. Via indium contacting, Au wires were 
attached as leads on samples at room temperature. Thus we were able to measure the 
transport properties of the fabricated devices by standard or quasi four- electrode 
method. Transport measurements were performed in PPMS (16T) and PHMF up to 
50 T. 
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